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The characteristics of the particle emitting source are deduced from low pt identified 
hadron spectra ((m-r — m ) < 1 GeV) and HBT radii using a hydrodynamic interpreta- 
tion. From the most peripheral to the most central data, the single particle spectra are fit 
simultaneously for all tt^, K^, and p/p using the parameterization in [ |TJ and assuming 
a linear transverse flow profile. Within the systematic uncertainties, the expansion pa- 



(N) 1 rameters Tf and Pt, respectively decrease and increase with the number of participants, 
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saturating for both at mid-centrality. The expansion using analytic calculations of the &t 
dependence of HBT radii in [ EI is fit to the data but no \ 2 minimum is found. 



g ■ 1. INTRODUCTION 

Identified charged hadrons in 11 different centrality selections [ [| and the transverse 
momentum (jpx) dependence of HBT radii in 9 k? bins [|] are measured in Au-Au collisions 
at 200 GeV by the PHENIX Experiment [ §. In both the 200 GeV and previously 
measured 130 GeV data [ |6|, the (j>t) of all particles increases from the most peripheral 
to the most central events and with heavier particle mass (mo). The dependence of the 
(Pt) on m Q suggests a radial expansion, and its dependence on the number of participant 
nucleons (N part ) may be due to an increasing radial expansion from peripheral to central 
events. The kx dependence of the HBT radii was also observed and interpreted as a radial 
expansion. 

Both the spectra and the kx dependence of the HBT radii are fit using parameteri- 
zations based on a simple model for the source, where fluid elements are each in local 
thermal equilibrium and move in space-time with a hydrodynamic expansion [[EH]- The 
assumptions are: (1) no temperature gradients, (2) longitudinal boost invariance along 
the collision axis z, (3) infinite extent in space-time rapidity rj, and (4) cylindrical sym- 
metry with radius r. The particles are emitted along a hyperbola of constant proper time 
tq = \ft 2 — z 2 and short emission duration, At < 1 fm/c. 

The mT dependence of the yield m ^ tdy \ y =o is calculated after integrating the source 
over space-time (azimuthal and rapidity coordinates) [ HJ . It is assumed that all particles 
decouple kinematically on the freeze-out hypersurface at the same freeze-out temperature 
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Figure 1. Simultaneous fits in the range Figure 2. The expansion in each centrality. 

(m T -m ) < 1 (solid lines) for tt~ (top) The top panel is T fo and the bottom is (3 T , 

and p (bottom) in all 11 centralities both plotted as a function of N part . 
(scaled for visual clarity) [ ||. The 7r 



resonance region is excluded in the fit. 

Tf , and that the particles collectively expand with a velocity profile Pt{^) = Prr/R 
where R is the geometric radius, r is the transverse coordinate, and (3t is the surface 
velocity. (For a box profile, the average velocity is ((3?) = 2/3^/3 [0). The particle 
density distribution is assumed to be independent of the radial position in the fits to the 
single particle spectra. In the previous 130 GeV analysis [f|, a Gaussian density profile 
increases (3t by ~ 2% with a negligible difference in Tf a , while a parabolic velocity profile 
increases (3t by 13% and Tf Q by 5%. 

We use analytic expressions to calculate the HBT radii [ |2|]. A linear flow rapidity 
profile in the transverse plane is assumed and a Gaussian distribution is used for the 
particle density dependence on r. The parameters are the geometric radius R, the freeze- 
out temperature T, the flow rapidity at the surface t]t (Pt = tanh (??t)) and the freeze-out 
proper time r . 

2. RESULTS 

2.1. Fitting the single particle spectra 

In order to minimize contributions from hard processes, all rnj- dependent particle yields 
are fit in the range (m^ — mo) < 1 GeV. As resonance decays are known to produce pions 
at low pt [ we place a lower pt threshold of 500 MeV/c on 7r in the fit. A similar 
approach was followed by NA44, E814, and other experiments at lower energies. In Fig. [I], 
7r~ and p yields are shown as a function of pr for each event centrality [ [|. The top 5 
centralities are scaled for visual clarity. The solid lines are the simultaneous fits in the 
limited px range. Similar results are obtained for K ± , tt + , and p yields. 
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Figure 3. Constraining fits to 
27t~ HBT radii in 10% cen- 
tral events using the expan- 
sion measured from the spec- 
tra. The shaded region is the 
systematic uncertainty from (3t 
and Tf . 



The systematic uncertainty in Tf a is determined by adding in quadrature the change 
in inverse slope due to the pt dependent uncertainties in each particle yield at low pr- 
For ti ± , K ± , and p/p, the uncertainty is ±10, ±13, and 16 MeV respectively. Added 
in quadrature, the total systematic uncertainty in the inverse slope is ±23 MeV. The 
systematic uncertainty in j3 T is dominated by the uncertainty in the p/p spectral shape 
at low pt and is determined by measuring the change in fi T after fitting for p T > 0.85 
GeV/c. The systematic uncertainty in j3 T is 17.5%. 

For the 5% most central events, particles are coupled to an expanding system with 
a surface velocity of (3t = 0.7 ± 0.2(syst.) and decouple at a common temperature of 
Tf = 110 ± 23(syst.) MeV with negligible statistical errors. For the most peripheral 
events, T fo = 135 ± 3(stat.) ± 23(syst.) and (3 T = 0.46 ± 0.02(stat.) ± 0.2(syst.). The 
statistical error only is included in the fit, resulting in x 2 /dof = 260.9/52 for the most 
central and 321.5/52 for the most peripheral events. At 130 GeV, similar results were 
obtained, with (3 T = 0.70 ± 0.01, T fo = 121 ± 4, and x 2 /dof = 34.0/40.0 for the most 
central events (statistical and systematic errors are added in quadrature before the fit) [ 
0. 

The fit results of all particles within each event centrality are shown in Fig. ^. The 
top panel is Tj Q and the bottom panel is Pt, both plotted as a function of N part . Within 
the systematic uncertainties, the expansion parameters respectively decrease and increase 
with the number of participants, saturating at mid-centrality. 

2.2. Fitting the k T dependence of the HBT radii 

The HBT radii are measured from identical charged n pairs in 9 kx selections for 10% 
central events [ |j. The systematic uncertainty in the data is 8.2%, 16.1%, and 8.3% for 
Rs, R , and i?L, respectively. A simultaneous fit to the HBT data could not be found over 
a broad range of parameter space. As an example, if the parameters (3t and Tf are set to 
the values from the spectra analysis, then the fit to the HBT results constrains R and Tq 
from the R s and i?L data respectively, yet the model overpredicts R Q by more than 3cr for 
all but the first uit data point (Fig. D). The systematic uncertainties in f3r and Tf Q are 
represented by the shaded region. Within these boundaries, R ranges between 6.9 — 16.8 
fm and r ranges between 11.2 — 16.7 fm/c. 

The x 2 contour levels of the expansion parameters Tf D (vertical) and rj^ (horizontal) 
are shown for simultaneous fits to the spectra and separate fits to each HBT radius in 
Fig. |j. We note that no x 2 minima are found, hence the contours are not closed. For the 
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Figure 4. The \ 2 contour lev- 
els for the expansion parame- 
ters Tf (vertical) and rjT (hor- 
izontal) after fitting 7r ± , K ± , 
p/p spectra and 27r HBT radii 
as indicated. Both scales are 
zero suppressed. 
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spectra, the contours are closed and show an anticorrelation, however there is no overlap 
with the HBT contours. The HBT radius R s prefers large flow rapidity tjt > 1.0 and low 
temperatures T fo < 50 MeV. The parameterization has the most difficulty reproducing 




3. CONCLUSION 

The single particle spectra are qualitatively described by a hydrodynamic parameteri- 
zation that assumes boost invariance and a linear transverse flow profile. The transverse 
expansion in 11 different centrality classes is extracted from the single particle spectra. 
Within the systematic uncertainties, the expansion parameters Tf Q and /?t, respectively 
decrease and increase with the number of participants, saturating at mid-centrality. Ex- 
pressions for the HBT radii based on similar hydrodynamic assumptions and Gaussian 
density profiles do not describe the identical ir pair data. Such fits worked well at CERN 
SPS energies [ HJ, but fail at RHIC energies. 
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